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The Yielding of the Earth to Disturbing Forces . 

By A. E. H. Love, F.R.S., Sedleian Professor of Natural Philosophy in 
the University of Oxford. 

(Received November 28, 1908,—Read January 14, 1909.) 

1. Any estimate of the rigidity of the Earth must be based partly on some* 
observation's from which a deformation of the Earth's surface can be inferred, 
and partly on some hypothesis as to the internal constitution of the Earth. 
The observations may be concerned with tides of long period, variations of 
the vertical, variations of latitude, and so on. The hypothesis must relate 
to the arrangement of the matter as regards density in different parts, and 
to the state of the parts in respect of solidity, compressibility, and so on. 
In the simplest hypothesis, the one on which Lord Kelvin's well-known 
estimate* was based, the Earth is treated as absolutely incompressible and 
of uniform density and rigidity. This hypothesis was adopted to simplify 
the problem, not because it is a true one. No matter is absolutely incom¬ 
pressible, and the Earth is not a body of uniform density. It cannot be* 
held to be probable that it is a body of uniform rigidity. But when any 
part of the hypothesis, e.g ., the assumption of uniform density, is discarded,, 
the estimate of rigidity is affected. Different estimates are obtained when 
different laws of density are assumed. Again, whatever hypothesis we adopt 
as regards the arrangement of the matter, so long as we consider the Earth 
to be absolutely incompressible and of uniform rigidity, different estimates 
of this rigidity are obtained by using observations of different phenomena. 
Variations of the vertical may give one value, variations of latitude a notably 
different value. It follows that “ the rigidity of the Earth " is not a definite 
physical constant. But there are two determinate constant numbers related 
to the methods that have been used for obtaining estimates of the rigidity 
of the Earth. One of these numbers specifies the amount by which the- 
surface of the Earth yields to forces of the type of the tide-generating 
attractions of the Sun and Moon. The other number specifies the amount 
by which the potential of the Earth is altered through the rearrangement of 
the matter within it when this matter is displaced by the deforming influence- 
of the Sun and Moon. If we adopt the ordinarily-accepted theory of the 
Figure of the Earth, the so-called theory of “ fluid equilibrium," and if we* 
make the very probable assumption that the physical constants of the 

* For Lord Kelvin’s work on the subject, reference may be made to his ‘Math. and 
Phys. Papers,’ vol. 3, or to Thomson and Tait’s ‘Nat. Phil.,’ Part 2, §§ 735—737 and 
833—847. 
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74 Prof. A. E. H. Love. [Nov. 28, 

matter within, the Earth, such as the density or the incompressibility, are 
nearly uniform over any spherical surface having its centre at the Earth’s 
centre, we can determine both these numbers without introducing any 
additional hypothesis as to the law of density or the state of the matter. 
We shall find, in fact, that observations of variations of latitude lead to a 
determination of the number related to the inequality of potential, and that, 
when this number is known, observations of variations of the vertical lead 
to a determination of the number related to the inequality of figure. 

[Note added, December 15, 1908.—This statement needs, perhaps, some 
additional qualification. It is assumed that, in calculating the two numbers 
from the two kinds of observations, we may adopt an equilibrium theory 
of the deformations produced in the Earth by the corresponding forces. 
If the constitution of the Earth is really such that an equilibrium theory 
of the effects produced in it by these forces is inadequate, we should expect 
a marked discordance of phase between the inequality of figure produced 
and the force producing it. Now Hecker’s observations, cited in §6 below, 
show that, in the case of the semidiurnal term in the variation of the 
vertical due to the lunar deflexion of gravity, the agreement of phase is 
close. If, however, an equilibrium theory is adequate, as it appears to be, 
for the semidiurnal corporeal tide, a similar theory must be adequate for the 
corporeal tides of long period and for the variations of latitude.] 

2. We take a to represent the mean radius of the Earth and p to represent 
the mean density. In the undisturbed state the Earth will be taken to 
rotate as if rigid about a principal axis. We shall denote by po,po, V 0 
density, pressure, and potential in the undisturbed state. In the theory of 
fluid equilibrium it is supposed that, for a first approximation, p 0 may be 
taken to be a function of r, the distance of a point from the centre. If po were 
known, V 0 would be known, and then would be determined by the equation 

1 dp 0 3V 0 __ 0 

po dr dr 

and the condition that po vanishes at r — a. In the same theory it is 
supposed that, for a second approximation, the surfaces of equal density, 
which are also equipotential surfaces and surfaces of equal pressure, may 
be regarded as oblate ellipsoids of revolution about the polar axis, the 
ellipticity E of the ellipsoids being a function of r, which is small for all 
values of r in the range a > r >0. Thus a surface of equal density, given, 
to a first approximation, by the equation r = K, is given, to a second 
approximation, by the equation 

r = {1+E(^—cos 2 #)}, 
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where 6 is the colatitude measured from the polar axis. Then the theory 
leads to the equation* 

( 1 ) 

where e denotes the value of E at r = a, co denotes the angular velocity of 
rotation, and g denotes the value of gravity at the surface, so that 

g = ^irypa, 

the constant of gravitation being denoted by 7 . Further, if C, A denote 
the moments of inertia of the rotating body about its axis of figure and 
a perpendicular axis passing through its centre, there is no difficulty in 
obtaining the equation 

°- A = TE\}4r^ d '- < 2) 

In calculating the constant of precession (C — A)/A, we may, with 
sufficient approximation, take A to be given by the formula 

3 A = 877 - f po r t A dr . (3) 

Jo 

The values of a, p, co, g, e, (C—A)/A are pretty accurately known by 
observation. 

This outline of the theory is sufficient for our purpose. We shall not 
need to take account of those further refinements in which quantities of the 
order e 2 are retained. 

3. Now consider the rotating body subjected to deforming forces. With 
a view to calculating the deformation we may disregard the effects of rotation 
and consider the undisturbed body as spherical and of density po at a poin 
distant r from its centre. Let the forces be derived from a potential which, 
in the region 0 < r <1 a, is expressible as a sum of spherical solid harmonics 
of positive integral degrees. Let W n be a single spherical harmonic term of 
this sum, the degree n of this term being a positive integer. Under the 
action of the forces derived from the potential W n the body will be deformed. 
Let U denote the radial component of the displacement, and A the cubical 
dilatation at any point P. The density at P in the disturbed state will not 
be the value of p 0 which belongs to the point P, but it will be 

po-U^-po&. (4 

Whatever theory of the constitution of the Earth may be adopted,! 
provided that the physical constants, such as density or incompressibility, are 

* Gf. J. H. Pratt, ‘ Figure of the Earth,’ 1865 edition, p. 79. 
t See the note added at the end of § 1. 
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Prof. A. E. H. Love. [Nov. 28, 


functions of r, the quantities U and A will be products of the spherical 
solid harmonic W n and some functions of r. The potential Y of the 
disturbed body at a point will not be Y 0 but it will be the sum of Y 0 and the 
potentials due to: (i) a volume distribution of density —U (dpo/dr)—poA 
throughout the region 0 < r < a, and (ii) a surface distribution of density 
[poU] r = a over the surface r = a. Hence the potential will differ from Y 0 by 
terms which are the products of the spherical solid harmonic W n and 
functions of r. In the most important cases n = 2, and we may write 

U = H(r)^, A=/(r)W? ; y = V a+ K(r) W* (5) 

9 9 

where H(r), /(r), K(r) are functions of r, and the constant 1 /g has been 
inserted in the expression for U for the sake of convenience. If, in fact, 
W 2 is a periodic term of the tide-generating potential, W 2 /g is the “ true 
equilibrium height ” of the corresponding tide, that is to say, it is the height 
of the harmonic inequality which the forces answering to W 2 would produce 
in an ocean covering a rigid spherical nucleus, of the same size and mass as 
the Earth, if the depth and density of the ocean were negligible. Erom the 
definition of K (r) we easily find the formula 


K(a) 


0 pcv 


-r«f(r)]dr. (6) 

In what follows we shall write h for H (a) and k for K (a), so that the 
equation of the disturbed surface is 

r = a + hW 2 /g , 

and the disturbance of potential at the surface is kW 2 . The numbers 
referred to above as specifying the inequality of figure and the inequality of 
potential are the numbers h and k. 

4. The horizontal pendulum is an instrument for measuring small 
displacements of the vertical. The deflexion is proportional to the 
gravitational attraction in a direction at right angles to the axis of the 
instrument. The attraction in question may be regarded as due to the 
Earth, the Sun, and the Moon. The axis of the instrument may be taken to 
be normal to the disturbed surface r = a +^W 2 /^. The forces acting on the 
pendulum are: (i) a force derived from the potential W 2 ; (ii) a force derived 
from the potential kW 2 ; (iii) the component of undisturbed gravity 
tangential to the surface r — a + hW 2 /g. If W 2 is the disturbing potential 
due to the Moon alone, the three forces act in the plane containing the 
Earth's centre, the Moon's centre, and the bob of the pendulum. We 
calculate the forces in this case, beginning with the force (iii). We have 

D 3 


W 2 = 


-r»(| cos*0'-i), 
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1908.] The Yielding of the Earth to Disturbing Forces. 77 


where M denotes the Moon’s mass, D the distance of the Earth’s centre from 
the Moon’s centre, r the distance of a point from the Earth’s centre, O' the 
angle between the radius vector drawn from the Earth’s centre to the point 
and the radius vector drawn from the Earth’s centre to the Moon’s centre. 
The surface r = a + hW 2 /g is approximately an ellipsoid of revolution, and 
its eccentricity e is given, correctly to the first order, by the equation 

2 _ 3toyM 


The angle between the central radius vector and the normal at any point 
of the ellipsoid is, to the first order, e 2 sin 0' cos O'. Hence the component of 
gravity tangential to the surface is, to the first order, 


2>ha^M. 

D 3 


sin 9' cos O', 


or 


h dW 2 
a SO' 


The sense of this force is that in which O' increases. The magnitudes of 
the other forces can be written down; and the resultant force acting at right 
angles to the axis of the instrument is 

l3W a kdW 2 hdW 2 

a dO' a dO' a dO' * ' 

If the Earth were absolutely rigid, the force would be a~ l dW 2 /dO'. Hence 
the ratio of the observed deflexion of the pendulum to that which would 
occur if the Earth were absolutely rigid is 


l + k-h:l. 

This result* has been obtained by regarding W 2 as the Moon’s tide¬ 
generating potential, but it is clear that the same expression would be 
arrived at if W 2 were any term of the tide-generating potential due to the 
action of both Sun and Moon, for the numbers h and k are the same for all 
disturbing potentials expressed by spherical solid harmonics of the second 
degree. 

5. An approximate expression for the height of any oceanic tide which 
follows the equilibrium law may be calculated by neglecting the depth and 
density of the ocean. Let W 2 be a term of the tide-generating potential. 
Then the corresponding inequality of the surface of the ocean, calculated in 
the way described, is (1+&)W 2 jg, and the inequality of the nucleus is 
hW 2 /g. Thus the height of the corresponding tidef is calculated as 


* An equivalent result for an incompressible sphere with a special law of density is 
given by W. Schweydar, £ Beitrage zur Geophysik,’ vol. 9, 1907, p. 41, but his argument 
is, in one place, so briefly expressed that I find it difficult to follow. 

+ This result is, of course, well known. 
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78 Prof. A. E. H. Love. [Nov. 28, 

-t-Jc—h)W 2 /g, and the ratio of the actual height to the true equilibrium 
height is 1+k — h : 1. Apart, therefore, from a correction depending on the 
depth of the ocean, the self-attraction of the waters, and the pressure of the 
tidal wave upon the nucleus, the quantity which is measured by means of 
observations of the long period tides is the same as that which is measured 
by means of the horizontal pendulum. It is, therefore, not surprising that 
estimates of the rigidity of the Earth which are based on observations of the 
fortnightly tide are nearly the same as those which are based on observations 
of the variations of the vertical. 

6. Lord Kelvin's estimate was based chiefly on observations of the fort¬ 
nightly tide in the Indian Ocean, from which observations it was inferred that 

l+Jc—h = f (8) 

very nearly. This result has been confirmed by W. Sehweydar,* by an 
analysis of much more numerous observations of fortnightly tides, and also 
by an analysis of several sets of observations made with horizontal pendulums. 
But the most striking confirmation is to be found in the investigations 
of 0. Hecker,f conducted by means of horizontal pendulums. His very 
precise and consistent results show that the equation (8) is remarkably exact. 

7. An analysis of many series of observations of variations of latitude 
led S. C. Chandler^ to the conclusion that such motions are roughly periodic 
with a period of about 427 days. If the Earth were an absolutely rigid body, 
having the figure of an oblate ellipsoid of revolution, with ellipticity equal 
to that found by geodetic observations, and having that ratio of moments of 
inertia about polar and equatorial axes which is deduced from the observed 
amount of precession, it would execute a free oscillation about the steady 
motion of rotation in a period of about 306 days, and the variations of 
latitude would have this period. The lengthening of the period from 
306 days to 427 days is due to the defect of rigidity. The result that a 
yielding of the Earth would lengthen the period was obtained by 
S. Newcomb,§ and the period was investigated theoretically by S. S. Hough.|| 
The theory has been placed in a very clear light by G. Herglotz.H 

8. We take the undisturbed axis of figure as axis of z, and choose axes of 
x and y in the equatorial plane. Let l, m denote the cosines of the angles 
which the instantaneous axis of the Earth's rotation makes with the axes of 

* Loc. cit. 

+ “ Beobachtungen an Horizontalpendeln ... /’ c Veroffentlichung d. k. Preuss. 
geodatischen Institutes/ No. 22, Berlin, 1907. 

X < Astronomical Journal/ vols. 11, 12, 15, 19, 21, 22 (1891—1902). 

§ ‘Mon. Not. K. Astr. Soe./ 1892. 

j) 4 Phil. Trans./ A, vol. 187, 1896. 

IT £ Zeitschr. f. Math. u. Phys./ vol. 52, 1905, p. 275. 
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x and y, The yielding of the Earth may be computed as the yielding of 
a sphere to forces derived from a potential W 2 , where 

W 2 = — co 2 (lxmy) z. (9) 

It is, therefore, expressed by a radial displacement U and a dilatation A, 
where 

U = H(r)W 2 /< 7 , A=f(r)W 2 /g. 

This displacement of the matter is to he superposed upon the inequality 
due to uniform rotation about the axis of z. The moments and products of 
inertia of the deformed rotating body with respect to the axes of x, y, z can 
be calculated. The angular velocity of the frame of moving axes of x, y , & 
has components la), m«, co referred to their instantaneous positions. The 
equations of motion of the body with the calculated moments and products 
of inertia referred to the moving axes can be formed by the ordinary methods 
of dynamics, and the period of free nutation can be deduced. This is 
Herglotz’ method, generalised by the inclusion of the effect of compressibility. 

9. It is found that, to the first order, the moments of inertia, A, C are 
unaltered, but that small products of inertia are introduced. The values of 
these products are given, correctly to the first order, by the equations 


\yzdm = 




\zxdm = 


8?r. 

15 
8t r 


m 


15 2 g 




dr. 


dr , 


( 10 ) 


j^cfan = 0, 


where the integrations in the left-hand members are taken through the mass 
of the disturbed body. By means of the result (2) we obtain the products of 
inertia in the forms* 


where 


-ma(C-A), -Z*( C-A), 

'd 


0 , 


« f?° Tr (, ' 5E) dr = Tg I} 0 \jr { " 6H W) -^ 6 /« 


dr. 


The period r of free nutation is given by the equation 

^ __ 27rA 

T “ ®(l-«) (C-A)' 


( 11 ) 

( 12 ) 

(13) 


* The formulae (11) and (13), and a formula equivalent to (12) in the case of incom¬ 
pressibility, were given by Herglotz (loc, cit.). The constant a is equivalent, in this case, 
to the constant which he denotes by v. The result expressed by (12) and (13) agrees 
with that obtained by Hough if the initial density p 0 is uniform and the matter 
incompressible. 
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Prof. A. E. H. Love. [Nov. 28, 


If the Earth were absolutely rigid the period t 0 would he 

2ttA 


To 


co (C—A) ’ 

so that the period is lengthened by the yielding in the ratio 1: (1 —a). 
Now if we introduce the results (1) and (6), we find 


(14) 


, To 7, « « 

'-r =“ = *^, 


€ — 


aco 4 

2g 


(15) 


It appears, therefore, that the constant k is determined. 

10. The period denoted by t 0 is known from the constant of precession 
to be about 306 days; the period denoted by t is known from observations 
of variations of latitude to be about 427 days. The value of a^jg is about 
2 ^-g-, and the value of e is about r • Hence we find 

k = t 4 5 nearly. (16) 

Since h—k = \ nearly, and k = T 4 y nearly, we find 

h = nearly. (17) 

These results may be expressed as follows:—The inequality produced in 
the potential of the Earth, near its surface, by the action of the Sun and Moon, 
is about four-fifteenths of the tide-generating potential, and the inequality 
produced in the surface of the Earth is about three-fifths of the true 
equilibrium height of the tide. The results hold for each of the partial 
tides answering to the several periodic terms of the tide-generating potential. 

11. If the matter within the Earth is assumed to be absolutely incom¬ 
pressible and of uniform density p, we should have k = |-A. If, further, it 
is assumed to be of uniform rigidity g, the theory of the deformation of an 
elastic sphere would give the result 

h = ( 18 ) 

\ gpaJ 

But the actual relation h—k = obtained from the tides and the varia¬ 
tions of the vertical, then gives h = or pu = -fV gpa. Since gpa = 3'5 x 10 12 
nearly in C.Gr.S. units, we have pi = 7*6 x 10 11 nearly. Thus the rigidity 
•calculated by this method is about the rigidity of steel (Kelvin’s estimate). 
If, instead of using the result h-k = l we use the value for h , we find 
from (18) the value gu = \gpa, or about 1*2 x 10 12 . The approximate 
agreement of this estimate of rigidity with that which has been deduced 
from the rate of transmission of the second phase of earthquake waves to 
great distances should not, I think, be regarded as of much importance. 

12. One way of deducing an estimate of rigidity from observations of 
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1908.] The Yielding of the Earth to Disturbing Forces. 81 

variations of latitude combined with the hypotheses of absolute incom¬ 
pressibility and uniform density and rigidity is to give to To, t, e, a» 2 / 2<7 in 
equation (15) their actual values, so that ^ ‘TF nearly. This method gives* 
h = nearly, and then, by equation (18), we find ya = 1*7 x 10 12 nearly. 
A second way is to give to t 0 /t and aa> 2 j2g their actual values, and to e the 
value 2 ^- 2 , which it would have, according to the theory of fluid equilibrium,, 
if the matter were homogeneous. This method* gives k = 0*4226 and 
fx = 9*2 x 10 11 nearly. A third way is to give to t and aco 2 /2g their actual 
values, and to t 0 and e the values which, according to the theory of fluid 
equilibrium and the theory of precession, they would have if the matter 
were homogeneous. We should then find &= 0*685 and = 4*4 x 10 11 * 
nearly. In view of these results the approximate agreement of Hough’s 
result with Kelvin’s estimate does not seem to have much importance^ 

13. In order to obtain an estimate of the rigidity of the Earth from 
observations of variations of latitude, it seems to be necessary to combine the 
observations with some hypothesis which will admit of the ellipticity e and 
the constant of precession (C — A)/A having their actual values. One of the 
simplest admissible constitutions is that which has been proposed by 
E. Wiechert.| He takes the Earth to consist of a solid nucleus of density 
8*206 enclosed in a solid shell of density 3*2, the ratio of the radius of the 
nucleus to the outer radius of the shell being 0*78. If we make the further 
assumptions that there is no slipping at the interface between nucleus and 
shell, and that the matter is absolutely incompressible and of uniform 
rigidity, we can calculate the numbers h and h in terms of the rigidity. If, 
then, the rigidity is adjusted so that the period of free nutation may be 
427 days, there results g, = 11*68 x 10 11 nearly. If, however, the rigidity is 
adjusted so that the displacement of a horizontal pendulum by tide¬ 
generating forces may be two-thirds of what it would be if the Earth were 
absolutely rigid, there results g, = 6*3 x 10 11 nearly. The first of these 
results was obtained by Herglotz by transforming the equations of elastic 
equilibrium for a body of variable density, and solving the transformed 
equations. The second result was obtained by Schweydar by using some 

* It is, in effect, the method adopted by Hough ( loc . cit.). 

+ This statement is not meant to suggest any doubt as to Kelvin's general conclusion 
that the Earth, as a whole, is a very rigid body. All the astronomical evidence confirms 
this conclusion. Kelvin's rather precise numerical estimate of rigidity, as some 7 or 
8 times 10 11 dynes per square centimetre, appears to be confirmed rather well by 
Hough’s estimate of about 9 times the same. Hough described his estimate as founded 
upon a “ reasonable hypothesis.” Two other reasonable hypotheses lead to and 17, 
where Hough’s leads to 9. These numbers confirm the general conclusion, not th& 
numerical estimate. 

X ‘ Gottingen Nachrichten,' 1897. 

VOL. LXXXII.—A. G 
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82 Prof. A. E. H. Love. [Nov. 28, 

analysis due to Chree.5 By an extension of the analysis which was applied 
to such problems by Kelvin, I have found the values 11*52 x 10 11 and 
6*3 x 10 11 . The discrepancy between the two estimates of rigidity, deduced 
respectively from the period of free nutation and the displacement of 
horizontal pendulums, shows, as we should expect, that the assumption of 
nniform rigidity is untenable. Schweydar has proposed to assume that 
there is one uniform rigidity for the nucleus and another for the shell, and 
finds that the results of the two kinds of observations can be reconciled by 
taking— 

fi for the nucleus = 2*02 x 10 12 , 

fji for the shell = 0*9 x 10 11 . 

The values which I have found are 2*15 x 10 12 and 0*86 x 10 11 . The slight 
disagreement between my results and those of Herglotz and Schweydar is 
due to slight differences in the assumed data. The general conclusion, that 
the rigidity of the nucleus may be much greater than that of ordinary 
materials at the surface, and the rigidity of the shell smaller than that of 
most rocks,f is more important than the numerical values. 

14. The results which have just been described may be taken as indicating 
the effect of heterogeneity on the estimate of rigidity. It appears that 
increase of density towards the centre compensates to some extent for defect 
of rigidity, and that increase of rigidity towards the centre can compensate 
for a considerable defect of rigidity in the superficial portions. The effect 
of compressibility is not known,]: but it seems improbable that the yielding 
of a compressible sphere with an assigned rigidity should be less than that of 
an incompressible one. The result that the rigidity of Wiechert’s shell may 
be less than that of most surface rocks led Schweydar to adopt Wiechert’s 
suggestion that there may exist a plastic sheet between the nucleus and the 
shell. I think it may be regarded as certain that there is not within 
a depth of 1400 km. a continuous layer of molten matter, separating 
the inner portions of the Earth’s body from the outer portions, and behaving 
as a fluid in respect of forces of the type of tide-generating forces. In order 
that the astronomical motions may be performed as we know they are, and 
that the surface may not yield to such forces more than we know it does, 
the portions of the Earth which are outside such a sheet, if it exists, must be 

* £ Cambridge Phil. Soc. Trans.,’ vol. 14, 1889. 

t The rigidities of many kinds of granite and marble lie between 2xlO n and 3X10 11 
C.G.S. units. See the memoir by P. D. Adams and E. G. Coker, ‘Publications of the 
Carnegie Institution,’ No. 46. 

f My attempt to determine this effect in ‘Cambridge Phil. Soc. Trans.,’ vol. 15, 1900, 
does not now appear to me to be satisfactory. 
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much more rigid than we can reasonably conceive them to be. No amount 
of rigidity of the nucleus would enable us to satisfy the conditions. To 
illustrate this statement quantitatively it will be necessary to set down some 
of the analysis of the problem. 

15. We consider a sphere of gravitating matter stratified in concentric 
spherical layers. The matter in any layer is taken to be of uniform density 
po and rigidity y, these quantities varying from layer to layer, and it is 
taken to be absolutely incompressible. If we discard the assumption of 
incompressibility, the problem becomes much more difficult. We suppose 
the sphere to be deformed by body forces derived from a potential W n , which 
is a spherical solid harmonic of the nVa degree. The initial pressure and 
potential Y 0 in the layer become po+p and Y 0 + Y n + W n , where Y n is the 
potential due to the inequalities at the disturbed boundaries of all the layers. 
Let (u, v, w) denote the displacement at any point (x, y, z) within a layer. 
The equations of equilibrium are three equations of the type 



'3Yo + aw + aw \ 



(19) 


dx dx ox 


and, since the undisturbed state is one of equilibrium, we can replace these 
by equations of the type 


^£+^ u + po ^ ( y n+Wn) = o. 


( 20 ) 


We have also the condition of incompressibility 

— o 


( 21 ) 


dx~*~dy + dz 

To solve these equations, we put 



( 22 ) 


where i is an integer, positive or negative, v and w are to be derived from u 
by replacing m successively by and Wi, and x successively by y and z ; also 
vji are spherical solid harmonics of the ith degree, and ^ is a spherical 
solid harmonic of the i\h degree defined by the equation 



(23) 


The summation refers to all the values of i that need be taken. 

16. With a view to the expression of the radial displacement U and the 
components of traction across a spherical surface r = E, it is convenient to 
introduce a spherical solid harmonic of the ( — i— 2)th degree defined by the 
equation 



(24) 
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Then we find rTJ = — 2 


2i-l 

L2(2i+1) 




/p2i+S 

2*+l 


0-t-aJ. 


(25) 


The value of U at r = E can be expressed by means of (25) as a sum of 
spherical surface harmonics of positive degrees, say SS j. We shall write 
for the corresponding sum of spherical solid harmonics of positive degrees, 
so that Q = ( r/T&y S j. If r = R is a boundary of a layer in the undisturbed 
state, and p 0 ', po" are the densities just inside and just outside this surface, 
the contribution of this boundary to the inequality of potential Y n is 




(26) 


according as the point at which Y n is estimated is within or without the 
surface r = R. The value of Y n is the sum of such contributions from all 
the bounding surfaces. 

The ^component of traction across the surface r = R is the value at r = E 
of the expression 


or 


- x - | 23o+/5o (Y„+W M )- A 62^if i | 


^•-i , r 
dx 


■ _ 9 ('hr 1 

2i +1 dx \r 2i ' 


2i +1 dx 


li+3 1 


</>- 


•<-*)}■ (27) 


The radial component of the traction is the value at r = E of the 
expression 


- {po+po (V.+w tt )} + pZ |(2 i +3) 


i 2 


+1 2i + 5 


^z+l 

2i 


2i + l 




(28) 


The tangential traction across the surface r = E is a vector quantity 
which has a definite magnitude and direction at each point of this surface, 
the direction being at right angles to the central radius vector. We may 
resolve this vector into components parallel to the axes. The component 
parallel to the axis of x is obtained by subtracting the ^-component of the 
radial traction from the ^-component of the resultant traction. We find 
that the ^-component of the tangential traction is the value at r = E of the 
expression 


£2 


t 




r 2 ctyz-l 


2^+1 dx 


r 2i + l 

2i-fl 0^\r 2i_1 / 


(2i+l) (2i + 3) 




2 i 


r 2i +3 


(2 1) (2i + 3) 


dx 


}■ 


(29) 
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17. The conditions to be satisfied at a disturbed boundary of a layer are 
conditions of continuity of displacement and traction. Except in the case of 
the normal traction, it is sufficient to form the conditions at the undisturbed 
boundary. Let letters with single accents denote the values of quantities 
such as Ui and fi just inside a boundary r = E, and letters with double 
accents denote the values of the corresponding quantities just outside the 
same boundary. The condition of continuity of the ^-component of dis¬ 
placement is 

R 2 ty ' <+1 , .R^'-i 

“2(iTT)-’- 1+ Tt & 


= u 


R 2 ty" i+1 


+ W -i -1 + 


R 2 ty'-i 


(30) 


2(^ + 1) dx 2i dx 

where we have picked out all the terms which contain spherical surface 
harmonics of the same degree i. We may now regard i as positive. The 
spherical solid harmonic of the ith degree, 

R 2 ( dyfr'j+i 
2(i+l)\ ox, 

1 r 2 * +1 fdf-j df'-A 


•u 


’ ) + ™Zi( w -t-1 —« -t-lj 


dx 


E 2i+ i 


+ : 


2 i E 2t 1 \ dx dx 
vanishes at the spherical surface r = E, and therefore vanishes for all finite 
values of r. We have two other like expressions which vanish in virtue 
of the conditions of continuity of the y- and ^-components of displacement. 
Hence, by the usual method, we deduce the two equations 


fi+i + 


v>214* 3 


i+X _5 

E*+ 5 2 E 2i+3 


and 

y£i + 3 


f-i -2 = same with doubly accented 
letters, (31) 


x* _2i+1 1 , 

E 2 9 -i-2 2 E 2i+ 5 + 


(2i + 3) 2 r 2i+ * 
4 E 2l+3 


f -i - 2 = same with doubly accented 
letters. (32) 


The condition of continuity of the ^-component of tangential traction is 






-U -i-\- 


^+3 


i.2i+l 


E^+i 

3 /jr'i-A i + 2 


2i+5 

i a 


» +1 


-2 r 


2i+l 


i-1 

2i +1' & V r 2 * -1 / 27+1 R K+1 a® 


(2t+l)(2*-l)R* +i aB 




0a; 2^ — 3 R 2 ‘~ 

(jJK+ty'-,-,) 

3 0 

(2i+l)(2i + 3)0a; 

2 0-2) r* +1 

1) R 2 0£c 

2 Q+3) 


Bi/r'-i 

dx 


(fV+Y-i.i) 


(2i-3)(2i 

+ 


(2t + 5) (2i+3) 


R2i+3 


5«j>'. 


dx 
(33) 


same with doubly accented letters. 
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Prom the three equations (two independent) of this type we obtain in the 
usual way the equation 


% (Jj -}- 9. W 9/7.4- K 3 


fx (i +1) (i + 3) ^i+i-T 



2i(i+-2)(2i + 5)r ai+8 
+ (2i + l)(2i + 3j“ E 2 


™2i+3 


2 (fr+1) (j+3) 


= same with doubly accented letters. (34) 


2i + 3 B 2i+5 J 


In forming the condition of continuity of normal traction at a disturbed 
boundary r — B + U, we may calculate all the terms that arise from the 
additional stress as if the boundary were the sphere r — E. The only terms 
in which we need take account of the displacement of the boundary are those 
that arise from the initial stress. The initial stress being hydrostatic 
pressure, the corresponding normal component of traction across any surface 
at a point is the pressure at that point with its sign changed. The initial 
pressure at a point on the surface r = B + U is expressible as the value of 
y? 0 + TJ (0y?o/9r) at a point on the surface r = E. When r = E is a boundary 
of a layer, p 0 is continuous in crossing the boundary, but dpojdr is not, for 
its value is po(dVo/dr), and p 0 is discontinuous. Hence the condition of 
continuity of normal traction at the surface r = B + U is 





> ) ^ 'i+\ 4- (2i + 3) 


2 i r 2i+s 2,(f+3) 1 

; i 1 "R2 * —i—2 O/i _L T?2i+ 5 


— same with doubly accented letters. 


2i +1 R? 9 * 2 2i + 5 E 2i+5 _ 


It is convenient to eliminate cf>'-i- 2 and (j>"-i- 2 by means of the equation 
(34). We thus obtain the equation 



+ 3 


{ 


2 i 2 + 8i -f- 9 

i.+ l 



^ f i+1 + (& -f 3) (2i-f ^'-i-2 



(35) 


= same with doubly accented letters. 


It is clear from these equations that in general the only functions of type 


<£, yjr which occur are </>„, yjr n , and the only functions of type m 

which occur are u n ~ h u n+ i, u~ n , u- n - 2 and the corresponding ^s and w’ s. In 
the central portion of the sphere the only functions which occur are those of 
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the types yfr n , cj)- n - 2 > 1 , u n + 1 . All these can be expressed in terms of V K 

by solving linear equations. It appears that U and Y n are multiples of the 
surface harmonic W n r ~ n . 

I obtained the results stated in § 13 above by applying the above analysis* 
to the case of a central spherical nucleus and an enclosing shell with 
Wiechert’s values for the densities and radii. 

18. We proceed to the example of an absolutely rigid nucleus, of density 
pi and radius a ly separated by a layer of fluid* of density p 2 from a solid shell 
of density p 3 , rigidity jn, and inner and outer radii a 2 and a 3 , the whole being 
subjected to body force derived from a potential W 2 , which is a spherical solid 
harmonic of the second degree. Let £", £ denote the spherical solid harmonics 
of the second degree to which the radial displacement U becomes equal at 
the surfaces r = a 2) r = a 3 respectively. Then at any point in the solid 
enclosing shell we have 



and 


The displacement at any point in the solid enclosing shell is expressed in 
terms of the functions (j> 2 , </>- 3 , fa, -\/r-. 3 . Equation (25) gives 



(38) 


The tangential traction vanishes at r = ct 2 and at r = a 3 , and hence we have 
the equations 



(39) 


The condition of continuity of normal traction at r = a 2 4- % is obtained 
from equation (35) in the form 



* In regard to the question of the adequacy of the statical theory to determine the 


behaviour of the supposed fluid layer, reference may be made to the Note at the end 


of § 1. 
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and the condition that the normal traction at r = a 3 4* £ vanishes is obtained 
in the same way in the form 

M*+ >+28 5'f-»- 8 $) 

~ 4 -7Pa ™ i (p 2 - ps) ~ T+ ipM -|1 (pi-P 2 ) 

+ (p*-/>»)^+paT««?l = 0. (41) 

19. It is sufficient to consider the case in which p 2 = ps, or the density of 
the fluid layer is the same as that of the solid enclosing shell. Writing g for 
the value of gravity at the surface and p for the mean density, so that 

g = ^TTypaz = 4 -tty -[(pi—ps)^ 


and putting 




W s 

9 


.af 1 
-Z + P2 

«3 S J 


^2 j ? 

a ?} ~~ 


(42) 


1 = h 


p gp2<* 3 


(43) 


we find that the value of h is given by the equation 

jll 456 -1800/ 3 +2688/ 5 -1800 )f + 45 6/ 10 “ 

24 + 40/ 3 —45/ 7 *— 19/ 10 
20. Suppose the fluid layer to be thin, and consider the case where the 
density is distributed according to Wiechert’s law. We have then p 2 = 3*2, 
p =s 5*58, / = 0*78. Remembering that gpa 2 = 3*5 x 10 12 nearly, and 
adjusting pu so that h may have its actual value f, we find fi = 3*5 x 10 i2 
nearly. It appears, therefore, that, even if the solid nucleus were absolutely 
rigid, and the enclosing shell were 1400 km. thick, the presence of a layer 
of fluid separating the nucleus from the enclosing shell w^ould increase very 
much the yielding of the surface. To prevent the surface from yielding 
more than it actually does, the rigidity of the enclosing shell would have to 
be nearly five times that of steel. If the enclosing shell were thinner, 
a still higher rigidity would be needed. For example, if it were 64 km. 
thick, and of density half the mean density, or about 2*8, the requisite 
rigidity of the enclosing shell, the nucleus being absolutely rigid and the 
fluid layer thin, would be about 50 x 10 12 dynes per square centimetre, or 
about sixty-six times the rigidity of steel. These numbers seem to me to be 
decisive against the hypothesis of the fluid layer. This conclusion does not 
negative the possible existence of areas of continental dimensions beneath 
which there may be molten matter. It means that such areas must be 
isolated; the molten matter beneath them cannot form a continuous sheet 
separating a central body from an enclosing crust. The conclusion does not 
negative the possible existence of a layer of comparatively small rigidity ; 
but, if there is such a layer, it must be rigid enough to prevent a finite 
slipping of the enclosing crust over the central body. 













